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There  is  given  a single  input-output  linear,  time-invariant 
plant  with  large  parameter  uncertainty  consisting  of  two  parallel 
branches,  one  of  which  has  n internal  sensing  points.  The 
objective  is  to  satisfy  specified  frequency  domain  bounds  on  the 
system  response  to  commands  and  disturbances  over  the  parameter 
range,  and  to  do  so  with  sensibly  minimum  net  effect  at  the  plant 
input,  of  the  n + 1 sensor  noise  sources.  The  basic  problem  is 
how  to  best  divide  the  feedback  burden  among  the  n + 1 available 
feedback  loops  L^.  The  procedure  developed  has  high  transparency. 
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giving  early  perspective  on  the  loop  bandwidths,  permitting 
approximate  loop  trade-offs  without  a detailed  design.  VJhile 
the  development  is  more  difficult  than  in  the  single  cascaded 
plant  system,  the  procedure  and  final  results  are  very  similar: 
Each  has  only  one  distinct  frequency  range  say  in  v/hich 

there  is  trade-off  between  and  and  ^ with 

steadily  increasing  loop  bandwidths  going  backwards  from  plant 
output  to  input.  It  is  shown  that  for  a class  of  problems  the 
sensor  noise  effects  can  be  tremendously  reduced,  when  compared 
to  an  optimum  single-loop  design  satisfying  the  same  specifica- 
tions. 
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a . 

1 

lower  bound  of  ki  (following  4) 

bounds  on  lT(ju))|  (1,  Fig.  2) 

b. 

1 

upper  bound  of  ki  (following  4). 

(oj) 

bounds  on  Li^(ja))  (Fig.  3,  II  Design  example) 

hf  boundary  on  Li^  (6b,  7,  Fig.  3) 

BW(L^) 

bandwidth  of  Li^  (IV) 

output  of  Pi  (Fig.  1) 

S 

system  output  due  to  D (Fig.  1) 

disturbance  input  (Fig.  1) 

system  function  (3,  8) 

®i 

excess  of  poles  over  zeros  (before  4) 

F 

prefilter  (Fig.  1) 

system  function  (4) 

compensation  function  (Fig.  1) 

hf 

high-frequency  (range)  (following  3) 

^i3 

j = A,B,C,;  i = 1, 2. . significant  u intervals 

^io 

i e,  hf  parameter  (4) 

4 

effective  hf  parameter  (4,  6) 

Li,  L° 

loop  transmission  (3,  6,  Fig.  4) 

^io 

nominal  value  of  Li  (II  Design  example.  Fig. 
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L.„  for  Xj  > (Fig.  3,  IV) 
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linear  state  variable  feedback  (I) 
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- <2,  Fig.  3) 

sensor  noise  input  (Fig.  1) 

as  sub  indicates  nominal  (Fig.  4) 

ith  plant  section  (Fig.  1) 

effective  plant  (3,  6,  12,  Al,  Fig.  1) 

system  command  input  (Fig.  1) 

as  sub  indicates  single  loop  (II,  Fig.  4) 

system  transfer  function  (1) 

variation  of  a set  (5,  9,  12) 

plant  input  (Fig.  1) 

trade-off  parameter  (Fig.  3,  IV) 

bound  on  M (2,  Fig.  3) 

= Pi/Pio  (9^  12) 

frequency  radians  per  sec. 

i = 1,2,...  significant  w values  of  (7,  Figs  3,  4) 

w at  which  Arg  = -ir  (Fig.  7,  Appendix  1) 


A SYNTHESIS  THEORY  FOR  A CLASS  OF  MULTIPLE-LOOP 


SYSTEMS  WITH  PLANT  UNCERTAINTY 
I.  INTRODUCTION 

There  are  two  distinct  approaches  to  the  design  of  multiple-loop 
linear  t ime- invar iant  systems.  One  of  these  LSVF  (linear  state- 
variable  feedback),  uses  the  optimal  quadratic  regulator  solution  and 
originally  secured  the  desired  system  poles  via  a constant  feedback 
gain  matrix,  driven  by  all  the  plant  states  [1,  2,  3 for  bibliography.] 
Later  refinements  were  observers  for  states  which  could  not  be  sensed 
and  of  prefilters  to  obtain  desired  zeros  l3,  ^] . The  problem  of 
parameter  uncertainty  is  being  currently  intensively  researched  using 
the  concept  of  "robustness"  [6-9].  LSVF  is  attractive  because  direct 
crank-turning  gives  a feedback  design  for  a multiple-loop  plant  of 
any  finite  complexity,  which  has  the  desired  poles  at  the  nominal  plant 
values  and  remains  stable  for  sufficiently  small  parameter  variations. 

A major  shortcoming  is  that  one  cannot  'design  to  specifications'  i.e., 
secure  specified  performance  bounds  over  a given  range  of  plant  parameter 
values.  Another  is  its  complete  neglect  of  the  price  paid  for  the  benefits 
of  feedback— the  bandwidths  of  the  loops.  Thus,  LSVF  insists  on  a feed- 
back structure  even  when  there  is  exact  knowledge  of  plant  parameters 
and  disturbances  - a situation  where  feedback  is  not  needed.  In  this 
sense  LSVF  is  a continuation  of  classical  network  synthesis,  using  a 
different  set  of  building  blocks,  because  its  primary  purpose  is 
pole-zero  realization  and  only  Incidentally  considers  the  uncertainty 
problem. 
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The  second  approach^ denoted  here  as  'quantitative  design'  is 
characterized  by  (1)  'design  to  specifications'  for  significant  plant 
uncertainty  and  disturbance  attenuation.  (2)  emphasis  on  loop  band- 
width minimization.  So  far,  these  have  been  secured  only  in  terms  of 
frequency  response  so  it  is  often  called  the  'classical'  approach, 
incorrectly  because  classical  control  theory  almost  completely  over- 
looked both  these  problems.  There  is  no  crank-turning  here,  but  purpose- 
ful design  for  sensitivity  reduction.  It  has  been  developed  only  for 

the  cascade  plant  structure  [10],  and  to  a certain  extent  for  the 

« 

multi-variable  two  matrix  deg ree-of- freedom  structure  [11].  This  need 
at  present  of  separate  development  for  different  structures,  compares 
unfavorably  with  LSVF  generality.  But  in  return  there  is  highly 
economical  design  to  specifications,  and  deep  understanding  of  the 
feedback  mechanism.  Also,  the  concept  of  ‘set  equivalence'  enables 
these  techniques  to  be  rigorously  applied  to  large  classes  of  linear 
and  nonlinear  uncertain  time-varying  systems  with  the  same  structures 
[i2,  13].  This  paper  extends  quantitative  design  to  the  cascade- 
parailei  mul tipie-ioop  structure  of  Fig.  1. 

Problem  Statement 

In  Fig.  1 the  P.  are  transfer  functions  of  sections  of  the  uncertain 
plant  and  N.  are  the  sensor  noise  sources  - drawn  heavy  to  emphasize 
they  are  constrained  and  unalterable.  There  is  independent  uncertainty 
of  the  parameters  of  each  P..  Despite  this  uncertainty,  the  system 
frequency  response  to  commands  T(ja))  = C ( jw) /R(ju))  is  to  satisfy 

specified  bounds 


0 < A^(tj)  < |T(jw)|  < 


(I) 
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It  has  been  shown  that  time-domain  bounds  on  the  output  and  its  deriv- 
atives of  any  order  [13]  can  be  achieved  by  satisfying  such  w-domain 
bounds.  The  problem  is  to  find  a sensibly  optimum  systematic  means  of 
dividing  the  'feedback  burden'  among  the  n + 1 available  loops.  All 
the  feedback  signals  go  to  the  plant  input  X , because  'plant  modification' 

[11]  is  assumed  not  allowed.  Thus,  in  Fig.  la,  X = C/P  and  each  C.  is 

determined  by  X and  the  P.  , so  the  C.  needed  to  obtain  a desired 

J ' 

art 

output  C independent  of  the  G.  . This  is  not  the  case  if  feedback  to 
internal  plant  variables  is  allowed. 

Disturbance  attenuation  is  another  major  reason  for  using  feed-  i 

back.  To  simplify  the  presentation,  for  the  present  the  only  require- 
ment in  Fig.  1 a,  (with  L defined  by  (3)  later),  is  that  |C/D[=((1+L)  ^1  < 
some  constant.  It  is  more  convenient  [15]  to  use 

M = Y a constant,  V u.  (2) 

Nonminimum-phase  plants,  unstable  plants  and  the  generality  of  the 
structure  in  Fig  la,  are  postponed  for  later  discussion,  except  to  note 
that  any  n + 2 degree  of  freedom  system  structure  [11]  may  be  used  e.g.. 

Fig.  lb  for  the  case  n = 2. 

The  above  is  a very  difficult  problem,  with  very  little  treatment 
In  general  and  none  at  all  for  Fig.  1.  It  obviously  does  not  lend  Itself 
to  a rigorous  mathematical  theorem-proving  treatment.  The  approach  taken 
Is  to  find  the  principal  design  factors  and  trade-offs,  based  on  the 
following  design  philosophy:  The  outer  loop  L from  C may  be  designed 
to  cope  only  with  the  uncertainty  in  P^,  P^,  which  can  give  an  L much 
more  economical  than  in  a single- loop  design  in  which  L must  cope 
with  all  P.  The  first  inner  loop  Lj  from  Cj  may  be  designed  to 
cope  only  v/ith  Pj,  with  possible  great  saving  compared  to  an  Lj  which 


h 


) 

r 

i 

t 
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copes  with  Similarly  the  second  inner  loop  need  cope  only 

with  P^,  etc.  The  result  is  considerable  transparency  and  insight, 
enabling  the  designer  to  decide  how  to  divide  the  feedback  burden  among 
the  loops. 

Simplifications  initially  made  in  order  to  concentrate  on  the 
essentials,  are  covered  in  Sec.  V. 

1 1 DESIGN  OF  OUTER  LOOP 


If  the  plant  is  minimum-phase  and  open-loop  stable  (1,2)  are 
achieveable  [15]  with  a single  loop  G.  = 0,  i = l,..,n,  G = 0 in 

Fig.  la.  But  the  resulting  = G^P  may  then  require  very  large  band- 
width, causing  great  amplification  of  the  sensor  noise  N,  as  in  a later 
example  - Fig.  5.  The  simplistic  approach,  later  justified,  is  to 
therefore  use  the  inner  loops  to  ease,  as  much  as  possible,  the  outer 
loop  burden.  In  Fig.  la,  let 


P = P P .•••PoP,,  P = P P.  + P 
a n n-1  21’  a b 


n 

P G. 


(1  + P G + P P ,G  , + P 

n n n n-I  n-1 

= ^ + Pg  ^^,(1+1),  L =1^^ 

C{s)  FGP  ^ L 

Rtil  " 1+(Gp7^)  " 1+L  • 

= L/P 


•PjGj)  + PG 


(3a-f) 


Fr  / 1 J.I  \ 


In  (3f)  the  sensor  noise  effect  is  examined  at  the  plant  input  X 
where  it  tends  to  be  large  [l6,  and  Fig  5 here],  causing  plant  saturation. 
In  the  high-frequency  range  (denoted  as  hf) , (3f)  L/P  where  1l(Jco)|«1 
but  I L/P I can  be  very  large  (Fig.  “ the  hf  range  is  the  major 
trouble  source.  Thus,  in  Fig.  5b  the  lowest  w range  with  large  and 
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sharp  peaking  of  |X/N|  is  300  rps  at  which,  from  Fig.  U,  the  nominal 
IlI/^  -l|6db.  Hence,  the  major  effort  in  sensor  noise  effect  minimiz- 
ation will  be  made  in  hf.  Since  P is  constrained  in  (4),  such  minimiz- 
ation requires  [l]  minimization.  But  from  (3)  L must  cope  with 
= P/^  uncertainty.  Therefore,  for  maximum  economy  of  L , choose 
the  G.,  i = i,...,  n in  of  (3c)  to  minimize  the  uncertainty  in 

P = P/^. 

^1^2”'^n^b‘^^c 

Consider  accordingly  the  uncertainty  in  P = 7^ p-  v - 

n n n 2 11 

e. 

hf  where  each  P.->-k./s  e.  the  excess  of  poles  over  zeros  of  P.. 

Ill  I 

Since  P,P„..P  P.  parallels  P in  Fig.  1,  it  is  assumed  that  (e,+e_+...e  ) 


I n 


1‘2'-  n b 

+ e.  = e +e,  = e . Hence,  at  hf 
babe 


_ k k. +k  . k 

pP^  abc  Ae  . 

P = — -> = — , where 

® s®<=(I+k  g +.,+k  g,]  s®c 

r,^r^  a 1 


n n 


. e A e +e  , 

-A  n - A n n-1 

G =gs  ,G  ,=g  ,s 
n ^n  ' n-1 


"n-r 


A ^a 
1 = 9,s 


(^a,b) 


The  range  of  k.  is  taken  as  [a.,b.],  b.>a.>0.  In  the  logarithmic 


complex  plane  (Nichols  chart),  P^  is  not  a point  but  a set  {P^} 

because  of  the  uncertainty.  For  any  fixed  kj,...,k^  values  the  set 

{P  } , due  to  {k.}  , {k  } in  (4a),  is  a vertical  line  whose  length, 
e DC 


Lgth  {P  }= 
e 


k b. +b 
a b c 


k a. +a 
a b c 


is  a function  of  k , and  is 
a 


db 


maximum  at  k = a If  b /a  > b. /a,  (at  b If  b /a  < b,/a,).  The 
3 3 cc*“bD  B cc““bb 

former  is  assumed  because  P^  is  in  parallel  with  “ see  Sec  V. 
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Hence,  due  to  a 1 1 the  k.  uncertainty  sets  kgth {P^ }>(a^b|^+b^)  / 

(a  a.+a  ),  with  equality  iff  3 g.  such  that  the  sets 

3 D O I 

{k  (k  ,k.  ,k  )}  of  (^»a)  c {k^(a  ,k  ,k  )}  as  the  k independently 

range  over  [a.,b.].  It  is  readily  seen  that  such  g.  exist  e.g. 

g_=...=g  =0,  b./b  sg  ^ 3. /a  , compatible  with  the  previous 
^ n Del  DC 

b /a  >b./a.  . (in  the  case  b /a  < b./a,  the  analagous,  compatible 
CC“bb  CC~DD 

condition  is  b./b  > g.  >a,/a  ).  Thus,  at  hf  the  best  the  inner 

D C “ I ^ D C 

loops  can  do  for  the  outer  loop  L,  leads  to  it  coping  with  a gain 
uncertainty  set 


Lgth 


a b. +b 
a b 


a d. +a 
a b 


£ 

c 


(5a, b) 


For  example  if  n=2,  all  a.=l,  b.  =40,  10,  60,  200  for  i = l,2,b,c 
then  in  a single-loop  design  must  handle  at  hf  {P}  of  length 

[(b  b.+b  )/(a  a.+a  )1j.  “Si-7  db  whereas  (5b)  given  1)2, 3db,  a saving  of 

dbe  dbCOD 

39.l*db, 


The  hf  region  is  most  important  for  sensor  noise,  and  the  hf  form 
of  P.  in  (4)  greatly  simplifies  the  problem  there.  But  design  of  the 
outer  loop  requires  the  uncertainty  set  for  the  entire  spectrum.  The 
complexity  of  the  calculations  for  general  P.  with  uncertain  poles 
and  zeros  would  obscure  the  important  features,  for  the  sake  (See.  V)  of 
a minor  point.  Therefore,  in  the  meantime  let  Pj  = k^/s"^  j = a,b,c 


so  (4,5)  apply  for  all 


Outer-loop  design  is  now  a single-loop 


problem  with  the  equivalent  plant  P^  of  (4a)  denoted  by 
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p P +P 

0 A ^aO  b 

1+P  G +...+P  ....P__P,_G, 

no  n no  20  10  I 


a k. +k 

a b c 

( l+a  g +. . .+a  g, 

' n^n  a I 


and  L°  = P°G 
e 


A 

A e 


(6a, b) 


The  super-oh  on  P^,  k^,  L indicates  L with  P^.P^.-.-P^  at 

their  nominal  values.  The  problem  is  to  find  G or  equivalently  a 

0 0 ^7 

nominal  = P^^G  so  that  (2)  and  Ailn |T( jo))  | < Jin  — 

are  satisfied.  The  optimum  design  for  this  single-loop  problem  [16]  is 

briefly  reviewed  here,  with  an  example  which  is  very  helpful  in 

explaining  the  multiple-loop  design  theory. 


Design  Example 

In  Fig.  la  let  n = 2 , P.  = k./s,  j = l,2,b;  P = k /s  ; 

J J c c 

aj  = 20,  a^  = 50,  a^  = I,  a^  = 1000,  bj  = 800,  b^  = 500,  b^  = 60, 

b^  = 200,000.  The  bounds  Aj  (co) , A2(w)  of  (1)  are  in  Fig.  2,  and 

Y = 2.3db  in  (2).  The  nominal  plant  values  are  taken  as  a.  (with  no 

loss  in  generality)  Note  1 : The  specifications  must  be  consistent 

with  physical  reality  i.e.,  it  is  crucial  [l6]  that3  w , such  that 

, A2(a))  0 

for  the  largest  variation  of  P^  <Jrn  • 'n  order  that  no 

sensitivity  reduction  be  needed  at  large  enough  u , permitting 
L(ja))  ->  0 as  0)  From  Fig.  2,  0)^“  55  for  a multiple  loop  design. 


i 

\ 

i 


? 


280  for  a single-loop  design. 
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Design  of  Outer  Loop 

The  procedure  in  [15]  is  follov^/ed.  From  (3a),  AJln|Tl  = in 


l + C 


L = PG=k/s^  of  (6)  with  k uncertainty  (5),  equal  to  ii2.3db. 
e e e 

Hence  in  the  Nichols  chart,  An  = An  + An  G is  any  point  on  line 

AB  in  Fig.  3.  of  length  i»2.3db.  By  means  of  G(ja))  this  vertical  line 

which  is  the  variation  set  can  be  translated,  but  not 

e 

rotated  to  any  region  in  the  Nichols  chart,  giving  the  variation  set  of 

L°.  V'{  L*^}  , with  the  nominal  ^ that  the  nominal 

L of  (3d),  denoted  by  and  the  nominal  P^^  of  (3d)  H.  P^^  , 

justifying  G = L^/P^q  = L°/pJ  . 

For  any  w say  w,  , one  finds  the  boundary  B(u,)  of  th^  set 

Lq(Jwj)  which  satisfy  (2)  and  An[T|  < An  example,  in 

Fig.  3 at  o)  = 20  , X,  is  satisfactory  for  L (J20)  because  the  range  of 

— -Q I is  from  A (M^  = -23.9  db)  to  B (M*^  = -.k  db)  giving  the 

A^{20) 

allowed  23-5  db  for  (Yo)  ' Similarly,  at  the  variation  is  from 

-22.7  to  .8  db.  Any  larger  |Lq|  at  the  same  Arg  is  satisfactory, 

but  not  smaller  |Lq[  . B(20)  is  thus  found.  Due  to  (1)  and  Note  1 

of  Design  Example,  as  o)  ->•  « B((i))  would  ->■  a vertical  line  at  -180“ 

extending  from  -A2.3  db  to  0 db.  But  (2)  gives  the  boundary  B of  Fig.  3 

H 


obtained  by  projecting  the  locus  of 


1+L, 


= Y=2.3  db  downward  by  ^»2.3  db. 


At  small  0)  , e.g. , w = .5  , (I)  dominates  - see  Fig.  3-  At  larger  u 

e.g.,  w = 2 , 10,  part  of  B(w)  is  due  to  (1)  and  part  is  due  to  (2). 

There  always  exists  w,.  such  that  for  to  > oi  , B(to)  = B„[l5,l6].  Here 

(o„  ' 70  rps.  The  generality  of  P at  hf  (i4a)  and  of  B^.  lead  to  a general 
H n 

shape  for  LQ(ja))  (Fig.  3,**)  in  large  hf  uncertainty  problems,  as  follows. 
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must  satisfy  B(a))  but  (3f)  at  hf  suggests  |LQ(ja))l  be 
decreased  as  rapidly  as  possible  vs  u.  As  s ->■  ^ k|^/s  . A 

reasonable  definition  of  optimum  is  one  satisfying  the  B(oj)  v/ith 

a minimum  for  a fixed  e^^.  Such  an  optimum  exists,  is  unique, 

lies  on  B((j)  at  each  u and  can  be  approximated  as  closely  as  des- 
ired by  a rational  function  [16].  There  is  trade-off  between  complexity 
of  the  rational  L_(s)  and  k,  . , so  a practical  sensibly  optimum 
LQ(ju)  is  as  shown  in  Figs.  The  shape  and  length  of  are 


important.  *-o  opt  tries  to  decrease  11^1  rapidly  vs  to,  but  in 

Fig.  3 B^  constrains  min  [Arg  LQ{jto)]  > - 130®,  with  corresponding 

minimum  average  d | L^Cjuj)  | /du)  = - (40)~  - 29db/decade  [ll].  Thus, 

ILqCJoo)!  must  decrease  rather  slowly  up  to  lo^  in  Figs.  3,*»  after 

which  the  permitted  decrease  of  Arg  Lp(jto),  at  bottom  of  B^^,  permits 

|Lp(jco)|  to  decrease  very  rapidly. 

This  paper  is  devoted  to  problems  where  3 an  [(0^,(0^]  interval 

I in  which  the  sensibly  optimum  LQ(jio)  has  the  shape  shown  in  Fig.  3* 

Plants  with  uncertain  highly  underdampled  pole-zero  pairs  (e.g.  bending 

modes)  could  be  included,  if  these  occur  at  lo  < o)^  and/or  to  > <o^. 

However,  the  multiple-loop  problem  is  complex  enough  without  bending 

modes,  so  this  class  is  omitted  here. 

It  is  seen  from  Figs.  3,^  that  the  hf  uncertainty  i.e.  of  k^  in 

(6),  is  the  factor  which  can  give  large  cost  of  feedback.  This  is 

because  the  length  of  B„  is  that  of  (k  }.  On  B , (l^I  must 

n C n U 

decrease  slowly  vs  to  while  [PqI  may  decrease  faster  and  at  hf  from 
(3f) , |X/n1  ^ * even  though  |Lq|<<  1 

(Figs.  A, 5).  Hence,  it  is  desirable  to  minimize  Lgth  (B^)  = Lgth  (k^). 
Use  of  an  inner  loop  permits  a maximum  reduction  of  39.^  db  here.  The 
saving  in  bandwidth  is  l»0/29  decades  ^0  f^'9*  ^)  • The 
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reduction  in  sensor  noise  effect  at  X is  enormous  (Fig  5b),  because  the 
rms  noise  value  is  obtained  by  integrating  arithmetic  values  on  a 
airthmetic  w scale. 

For  later  use,  the  following  o)  intervals  in  Fig.  3,^  are  emphasized: 

'a  " ~ [0,90),  Ig  ^ ~ ^90,  330) 

!(.  = ”)  ~ [330,  “>).  (7) 

The  design  of  the  first  inner  loop  Lj  is  decisively  influenced  by 
these  Intervals  of  the  outer  loop  L^. 


Ill  DESIGN  OF  INNER  LOOPS 

First  Inner  Loop  Lj 


In  II  the  inner  loops  were  apparently  sacrificed,  in  order  to 

obtain  the  most  economical  outer  loop  and  thereby  minimize  the  effect  of 

e 

sensor  noise  N at  X.  G.  = ...  = G = 0,  G,  = s b./b  were  found 

2 n ' 1 be 

satisfactory  for  this  purpose.  The  obvious  criticism  is  that  this  Gj, 
besides  being  impractical,  would  tremendously  amplify  hf  Nj  noise 
effect  at  X (Fig.  1)  and  likely  more  than  cancel  the  benefit  gained 
for  L.  The  answer  is  that  while  these  Gj  are  satisfactory,  there 
are  other  much  smaller  acceptable  values.  This  is  due  to  the  mechanics 
of  sensitivity  reduction  such  that  LQ(ju)j)  optimally  designed  to  handle 


an  uncertainty  set  V(a)j)  can  in  practice  handle  a set  Vy(h)j) 

much  larger  than  V(m,)  (e.g.  Fig  8 of  [10]).  So,  the  next  step  is  to 


find  the  bounds  Bj(a))  on  the  first  inner  nominal  loop  LjQ(jto)  such 
that  the  economical  of  1 1 is  satisfactory.  The  bounds  on  Ljq  are. 

In  fact,  very  modest. 


1 1 


For  this  purpose  (3c, f)  are  entended  as  follows.  Let 
Jb  =/^,(l+L)  = [(l+P^G^+...+P^...p2G2)  + P^G,](1+L)  = + P^G^Xl+L) 

Lj  = P^Gj 


A 


X 


L,/P 
I a 


N J^-  ^(T+L,){1+L)  (1+L,)(1+L) 


(8a-c) 


= L,/P  in  the  crucial  hf.  Hence  to  minimize  |X/N, I at  fixed  L , 
la  I * 

minimize  |Lj|  . But  Lj  must  cope  with  the  uncertainty  in  Pj  , P^.-.-P^ 
ignored  by  . However,  if  G^  can  cope  with  ^2’'‘''^n  *-]  need 


only  cope  with  P 


Lj  is  designed  accordingly  and  denoted  by  Lj  to 


indicate  its  neglect  of  p2,...,P^  uncertainty.  So  now,  P^  of  (7a)  is 
replaced  by 


1 A Pk+P 

pi  ^ I 20  no  b c 


X.P  P.P 
1 ao  b c 


e i+p^G+...+p  ..,p„-p,G,  “ ypTru^rcTj 

no  n no  20  1 1 «</ 20  1 10 


A ^ • 

' ^0 


(9a-c) 


with 


-iT'  -Vipt) . 


X.a  k,  + k 
lab  c 


‘ - .^iho 


Instead  of  of  (5a).  In  (9)  Lj^  is  the  nominal 

^?0  ~ ^10  ^0  “ '"O^  nominal  P.^,  for  i = 1,...,  n. 

Lq  was  designed  to  handle  P®  with  its  but  now  it  must  handle 

P^  with  its  What  are  the  bounds  Bj(aj)  o-  LjQ(ja))  so  that 

the  original  remains  satisfactory?  This  equestion  may  be  answered 

by  simple  trying  Ljq  values  and  checking  if  (1,2)  are  satisfied.  It 
is  found  that  the  Bj(u))  are  decisively  influenced  by  the  intervals 
*A’  *B’  ’c  *"0  results  are  stated  here  and  their 

explanation  in  Appendix  1. 
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Nature  of  Bj  (ui)  bounds  on  Ljq 

(1)  For  cael^=[0,90),B|(u))  arc  upper  bounds,  i .e.  , 

w)  [ must  be  < some  value  which  is  a function  of 
ArgL|Q(ju))  - Fig.  6a. 

(2)  For  u)  e l„  = [90,330),  B, (u)  are  lower  ones  precluding 

D I ■ ■ “ ■ 

Ljq  h 0 (Fig.  6b). 

(3)  For  0)  e Ij,  = [330,®),  Bj  (w)  in  Fig.  6b  are  closed  curves  in 

the  Nichols  Chart  which  tend  to  a vertical  line  of  length 


at  ArgL 


k ' db 


10 


u . 


Just  as  in  the  design  of  , so  the  optimum  would  lie  on 

B| (w)  at  all  u but  is  in  practice  approximated  by  a rational  function  - 
Figs.  1»,6.  One  may  define  intervals  of  similar  to  1^,  1^,  1^  of 


A’  ‘B’  C 


Lq  , i.e.,  in  Figs.  h.S:  I = [0,a.,-)  , 1,^  = [a,,-,  a.,  + ) , 1,^  = (o,+  *®) 

Here  l,„  = 0 because  B,.,  has  zero  width.  In  practice,  one  would 

Id  In 

likely  (in  addition  to  (2)),  assign  bounds  in  Fig.  1,  on 


i _ • 

D.  {1+L)tl+L^)...(l+L.) 


(10) 


leading  to  finite-width  B.^  and  larger  l.„  . Such  finite  E.,,  are 

In  Id  In 

easily  added  in  Figs.  6,7,  but  are  omitted  here  for  simplicity. 


Second  Inner  Loop 

The  above  discussion  is  repeated  for  L^t  but  now  uncertainty  is 

included  with  (8)  extended  to 


+ P^...p2G2l(l+L,)(l+L) 

- (^3  + P^...P2G2)(1+L,)(1+L)  =^(1+L2)(1+L,)(1+L), 
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An  2 2 


7^ 


-X  S 


L/P...P 
In  2 


^2  ^ (I+L2){I+L,)(I+L) 


(1 la-c) 


“ ^2^^n'‘’^2  crucial  hf.  To  minimize  the  latter  it  is  best  to  let 

handle  uncertainty  only,  leading  to  (cf  (9a)) 


•>  A nPu+P 

,2  A I 2 aO  b c 

e 1+P  G+...+P  ...P,  G,+P  ...P,„P^G,+P  .,.P,„P^P.G 

no  n no  3o  3 no  30  2 2 no  30  2 I I 


(A,X^P  P,_+P  ) 
' 1 2 ao  b c 


^ X = -^ 

30'^’‘‘'^2^20'‘'^1^2^10^''‘'^?0^^  ’ 2 P 


20 


(12a-c) 


wi  th 


{P^} 
e 


= < 


X,X„a  k,+k 
I 2 a b c 


instead  of  the  smaller  uncertainty  set  7/ ^ of  (9c). 

The  next  step  is  to  the  find  82(0)),  the  bounds  on  L^q  £ L2Q,  so 
that  L , L,-  designed  for  remain  satisfactory  for  'If^-  The 

• . r X 

resulting  B2(u)^  simi  lar  to  Bj(a)):  upper  bounds  in  liA=  lO.,o,^). 

lower  ones  In  l,„  and  closed  curves  merging  into  a of 

I B 


length  —1  - see  Fig.  7*  The  explanation  is  given  in  Appendix  2. 

\^Zl  db 

One  can  continue  indefinitely  in  this  manner.  The  resulting  L2Q 


(Fig.  1»,7)  has  three  intervals  12^.  I2B’  *20  decisively  influence 

the  bounds  on  a handle  P^  uncertainty,  etc.  The 

general  forms  for  the /5^  j,  L.,  etc.  are  for  i = l,...,n  ^ 


\k 

£\^  r ’ p - . . P . G . 

A -A*,*-,), 


„ L./P  ...P. 

-X  _ i n I 

N.  " d+U  )(1+L.  j)...(l+L) 


(13a-d) 


P = 


(A.X, 


.X.P  P.  + P ) 
I ao  b c 


Note  that  F in  Fig.  la  is  available  from  Eq.  3e  as  soon  as  is 

known,  by  associating  a nominal  with  the  nominal  ^q(s)  • But  G. 


is  not  known  until  L , L , 

n n-1 


_ ^no*^ n+1 ,0  _ . 


no 


/P  . 

no  no 


L.  are  known.  Thus  from  (1 3b) 


^n-l  ^n-1  .O'^nO^'^no^n-l  ,0 


with  A/  = I + L etc. 
no  nO 

Generality  of  structure.  In  the  system  considered,  input  R in  Fig.  1 

and  n+1  plant  outputs  are  available  for  processing,  permitting  an 

infinitude  of  n + 2 degree  of  freedom  structures  [11].  The  n + 2 

fundamental  system  functions  are  the  system  transfer  function  T(s)  = C/R 

and  the  n + 1 loops  L,  Lj,...,L^.  In  any  acceptable  structure,  L Is 

gotten  by  cutting  the  outer  loop  just  after  the  C sensor,  giving  in 

Fig.  lb,  L = P(iHHjH24(?i , *=  1 + P2H2  + first 

cut  and  with  another  cut  after  the  Cj  sensor,  gives  Lj  *=  P2^|^1^2  , 

“ 1 + P2^2"  always  of  the  form  T = ij;L/(l+L),  i|;  independent 

of  P.,  4’  = l/H  In  Fig.  lb.  The  design  technique  provides  T and 

the  nominal  L from  which  the  compensations  G.  (of  Fig.  la)  or 
iO  • 

H.  (of  Fig.  lb)  or  those  of  any  other  structure  are  derived.  The 
excess  of  poles  over  zeros  assigned  to  T(s)  e^^  must  be  compatible 
with  the  structure.  In  Fig.  la,  e^  *=  e^.  + e^^  **  e^  + Cp  + e^  each  a ^ 

positive  integer  but  in  Fig.  lb,  e^^  « Tc.  «=  Cj  + e^,  i = Q,H|,H2,P  H. 
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IV  PRACTICAL  DESIGN  PROCEDURE  AND  TRADE  OFFS 


Sections  II,  III  described  a design  procedure  based  on  the  best 

(most  economical)  , subject  to  the  best  •••  » subject  to  the 

best  Lj  , in  turn  to  the  best  L ; first  preference  is  given  to  L , 

then  Lj  , etc.  This  section  shows  how  II,  III  provide  the  perspective 

for  making  reasonable  trade-offs  between  the  loops  early  in  the  game, 

without  a detailed  design.  The  display  in  Fig.  h is  used.  The  first 

step  is  an  approximate  single-loop  L^.  design.  The  low  frequency  bounds 

Bg(ti))  based  on  P of  {3b)  are  used  which  are  hardly  different  from 

those  based  on  P^  of  (7a)  - see  Sec.  V.  There  is  no  need  for  a detailed 
e 

design  of  L-_  for  u > u , at  which  reaches  B,,-  (analogs  of 

dU  QS  oU  ho 

Wj  » B,,  in  Fig.  3)-  Thus,  the  slope  of  L-„  is  known  on  and  the 

d n oU  no 

length  of  B„c  is  that  of  the  hf  uncertainty  of  P . The  slope  of  L.,. 

Ho  so 

for  w > 0)  - (analog  of  w ) is  the  same  as  of  L-  for  ui  > ui  -fcf 

Lq  t Lgp  in  Fig.  A).  Having  , the  approximate  is  immediately 

available  because  B^  is  known  {39-^t  db  shorter  than  ^^15)  • 

Next,  sketch  an  approximate  L,.  as  follows.  Il.^I  is  near 

0)^  and  its  approximately  value  is  obtained  by  the  method  of  Appendix  1, 

Fig.  Alb.  The  shape  of  |l-|ol  ^ fairly  standard.  Its 

slope  is  « -30  db/decade  from  ui^  to  in  Fig.  If  until 

|LjQ(ju)|  = 20  log  - A db  is  attained  (A  a small  gain  margin),  after 

which  it  is  ~ constant  for  1-1.5  octaves  ((o,  to  Q.)  , followed  by  a 

I X ^ 

slope  of  -20  db/decade  , with  e^^^  the  chosen  excess  of  Lj^  poles 
over  zeros.  Analogous  to  , L*q  (Fig.  *>)  is  coping  with  all  of 

P uncertainty.  Similarly,  an  approximate  is  obtained.  Il.-I 

is  between  w,  and  Q,  (Fig.  k)  and  its  value  can  be  found  from  Appendix  2. 

I X j 

For  u > Qj  , Its  shape  is  similar  to  that  of  |1-^q1  for  to  > . For 


1 
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n > 2 , the  procedure  is  continued  with  Il.^I  near  , etc. 

’ ' 30 'max  2x 

The  next  step  is  to  sketch  (Fig.  k)  |P^|,  IPjo'-'^nol’  l^20'’'^nol’ 

The  noise  amplifications  in  hf  (3f»  13)  are  I^-jq/Pq!  » ’ 

respectively,  easily  obtained  by  subtraction  of  the  db  values.  The  sensor 
noise  effect  |X.|  is  gotten  by  multiplying  |X/Nj|  hy  |n.|  . 

Trade-offs  between  the  L.  are  now  considered,  e.g.,  vs.  Lj^  . 

Lq  of  I I is  one  extreme,  is  the  other  and  intermediate  designs  are 

possible.  One  poorer  by  a = 5 db  is  shown  in  Fig.  3.  postponing  the 
'^l'^2'^3  P'9-  ^ until  |Lq|  is  less  by  5 more  db  with 

> (‘‘>2^01  d*  return,  the  peak  of  the  new  I*-jq|  (Appendix  1,  \ 

Fig.  Alb)  is  ~ -18.5  instead  of  -9db.  Trade-off  between  and  Lj 

is  made  with  no  reference  to  L^,  L^,....  Trade-off  between  Lj^  and 
L^q  is  done  in  the  same  manner  etc. 

Bandwidth  Propagation  and  Similarity  with  the  Cascade  Plant  structure 

Let  the  bandwidth  BW(L.)  be  arbitrarily  defined  as  that  at  which 
Il.qI  achieves  its  final  asymptotic  slope;  for  L^,  for  Ljq, 

for  L2Q  in  Fig.  k.  BW(L.)  increases  with  i.  This  phenomenon 
occurs  in  precisely  the  same  manner  in  the  cascade-system  llO].  The 
relations  between  the  L.  , the  role  of  b./a.,  the  sensor  noise  effects 
and  trade-offs  etc.  are  very  similar  in  the  two  structures.  However,  the 
values  of  Il.-I  , are  different  and  the  derivation  is  more  difficult 
here.  Here,  at  each  new  L.  stage,  one  must  use  a more  complex  form 
of  P^.  In  the  cascade  system  the  step  from  i to  i -F  1 is  identical 
the  that  from  i - 1 to  i.  But  the  final  results  are  remarkably  similar. 

In  Fig.  ^2  comparable  with  BW(LgQ)  at  X^,  a 

little  larger  due  to  the  extra  few  db  of  gain  margin  needed  per  section. 


mamtmmmm 
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Thus,  the  final  cut-off  frequency  for  a single-loop  design  is  comparable 
to  that  for  a multiple-loop  design,  but  they  are  associated  with  different 
loops  so  there  can  be  a great  improvement  in  sensor  noise  effect.  Thus, 
in  Fig.  = -22  + 91  = 69  db,  while  |P2oldr  l^oldb  " 

an  improvement  if  ^ 127“69  = 58  db.  In  practice  it  is  reasonable 

to  assume  that  the  plant  power  levels  and  with  them  the  sensor  noise 
levels  Increase  in  preceding  from  input  to  output.  The  design  procedure 
is  highly  transparent  permitting  a good  estimate  of  the  optimum  division 
between  the  feedback  loops,  without  a detailed  design. 


High-frequency  uncertainty 

Clearly,  multiple-loop  design  can  be  highly  superior  to  single-loop, 


for  large  hf  plant  uncertainty.  The  linearized  plant  model  is  usually 


due  to  linearization  of  a nonlinear  about  an  operating  point  or  trajectory. 

K 


Large  variations  can  exist  due  to  different  operating  points,  e.g.  in 
flight  control  [17],  where  values  > 1000  have  been  reported. 

It  has  been  proven  that  in  a large  class  of  linear  and  nonlinear 
time-varying  uncertain  plants  the  latter  can  be  represented  for  synthesis 


purposes  by  an  equivalent  linear  time-invariant  uncertain  plant  set  P [s] 

eq 


[12,  13].  The  set  equivalence  is  exact  with  respect  to  a prescribed 
acceptable  plant  output  set.  Linear  time  invariant  design  applied  to 


the  problem  is  guaranteed  to  work  for  the  original  nonlinear 

W/  /JO 

problem.  A nonlinear  plant^can  thus  generate  large  hf  uncertainty  in 


P^q[s],  e.g.  consider  y = k^x^  , x the  input  and  y the  output. 


Suppose  fairly  linear  response  is  desired  for  y = A^(l-e~*)^,  A e [0.5,5] 


To  find  P [s] , evaluate  Y(s)  „ 6kA'^  ....  ».  / 

eq  W ° ^ •(s+2)'(s-k3)  '"‘disease.  Smee  I 


A e[0.5,5],  the  hf  gain  of  P varies  by  a factor  of  100,  due  to  A^. 

eq 


For  a simple  dynamic  example,  consider  ^ + By*^^  sgn  y = kx,  giving 


14 


p = — 

(s+3) [BAs+6A^+2BA] 
of  100. 
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6kA^  2 2 

->•  6kA  /Bs  at  hf,  with  uncertainty  factor 


V.  JUSTIFICATION  OF  ASSUMPTIONS 

General  plants.  This  section  is  devoted  to  the  justification  of 

e.  \ 

simplifying  assumptions  in  II,  III.  One  was  use  of  P.  = k./s  for  \ 

' ' \ 

all  0),  not  just  in  hf  where  it  is  applicable.  Recall  in  Sec.  II  the  \ 


\ 


first  step  was  to  find  the  smallest  {P  } of  (6a),  by  minimizing  over  / 

® / 
G ,...,G  and  the  values  of  P,_,...,P  Suppose  P.  = k./(s+q.) 


/ 


...V,  V,.  .jQ,...,.^Q.  ..  _ j. 

with  k j uncertain.  This  minimization  problem  is  extremely 

difficult  at  medium  o).  Fortunately  it  makes  little  difference  if  it 

is  not  done  at  all.  The  reason  is  that  which  made  unnecessary  in 

I . I in  I,,  etc.,  i.e.  under  certain  conditions  there  is  little 

A’  20  . lA 


difference  in  L-=GP 


. needed,  whether  {P}  = set  S,  or  set  S-«S,  . 
0 o'min.  1 21 

In  Fig.  3,  suppose  that  instead  of  AB  (A  at  X^) , the  uncertainty  set 
Is  ABEFG  v/ith  E,F  extending  even  to  at  results  in  almost  the 

same  AJlnlT|  for  both  (23.85  db  instead  of  23.5db). 

It  is  therefore  concluded  that  in  most  of  1^,  {P}  of  (3b)  be  used 
for  Lq  design,  just  as  In  design.  P^  is  used  only  for  u where 

Pj  is  well  approximated  by  kj/s  •*.  This  has  been  verified  for  several 
numerical  examples;  e.g.  for  n = 1 with  P^  = k/(s+q^)^P^  = 

P = k /s(s+q  ),  k^  e[l,l»00]  k.  e[l,60],  k^  e[l,200],  q e[0.5,2], 
all  independently  uncertain.  The  maximum  difference  in  the  two  B(a)) 
is  Only  three  db  even  though  the  difference  between  {P}  and  {P®}  is  ~ 
i*0db.  If  this  conclusion  is  incorrect  for  an  unusual  case,  then  it  is 
also  likely  that  the  obligations  on  Ljq  in  1^  will  be  greater  too 


By  using  (P)  in  medium  u,  one  is  certain  that  the  obligations  on 


cJ 


'I 


f 
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in  1^  will  be  negligible,  as  in  Sec.  II.  The  simple  and  transparent 
forecasting  of  Sec.  IV  may  then  be  used.  If  these  indicate  less  than  de- 
sired saving  in  sensor  noise  effect,  then  one  can  return  to  check  if 
greater  saving  is  possible  with  in  1^  . 

Another  assumption  in  II  was  b /a  > b /a.  . If  the  opposite  is 

c c b b 

true  then  minimum  Lgth  {P^}  is  at  k = b of  value  (b  b,+b  )/{b  a, +a  ) . 

^e  aa  abcabc 

There  exist  a set  of  g.  which  achieve  this  and  the  procedure  is  precisely 
the  same  as  before.  A third  assumption  is  that  e + e.  = e giving  (^) 

d D C 

with  {P  1 in  hf  a vertical  line  in  the  Nichols  chart.  If 
e 

|e  + e,  - e 1 = 6 5^  0 is  even,  the  result  is  also  a vertical  line 
'a  b c 

whose  length  is  a function  of  u.  The  design  procedure  is  basically  . 


the  same.  It  is  possible  that  P{jw)  = 0 at  finite  w = at  some 
combinations  of  parameters,  giving  T(ja)j)  = 0.  If  so,  the  specifica- 
tions on  T(jto)  and  C/D(jai)  must  allow  for  this.  If  6 is  odd, 
design  is  more  complicated  because  min.  Lgth  {P^}  does  not  necessarily 

exist.  The  range  of  {P^}  (i.e.,  at  any  fixed  P value)  is  no 

G d 3 

longer  a line  but  a two-dimensional  region  and  there  may  not  exist  a 

set  of  g.  values  in  (4a)  such  that  the  resul t i ng  1/{P^}  fits  into 

• Q e a 

any  one  It  is  then  a matter  of  judgment  how  to  exploit  the 

available  freedom  to  optimize  Lq.  This  case  has  not  been  studied  in 
detail.  However,  the  design  technique  of  secs.  II,  III  provides  the 
understanding  for  good  use  of  the  design  variables.  One  knows  the  kinds 
of  distortions  of  the  uncertainty  set  which  are  useful  in  relation  to 

Another  assumption  was  that  the  disturbance  attenuation  was  a minor 
problem,  dealt  with  by  (2,10).  The  procedure  is  basically  the  same  if 
it  Is  a major  problem,  for  then  = C/D  must  satisfy  | Cp( jto)  | <Y(aj)  over  j 
{P}.  This  can  be  translated  [I5l  into  bounds  ^^(w)  on  LQ(ju)).  The 
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more  stringent  of  Bp(w)  and  of  B(w)  due  to  (1),  is  used  but  there- 
after the  design  procedure  is  the  same. 


Unstable  and  Nonminimum-phase  plants.  Open-loop  stable  minimum-phase 

plants  were  assumed  in  II,  III  for  simplicity.  But  clearly  the  design 

procedure  applies  so  long  as  the  L.q  exist  which  satisfy  the  B.(a)/. 

Consider  first.  It  must  handle  giving  (Sec.  II)  a single-loop 

problem.  The  latter  is  solvable  if  (Pg)  contains  open-loop  poles  whose 

range  of  uncertainty  includes  part  of  the  right-half  as  v.(e!  1 as  the 

left-half  plane  [15,18].  If  however,  {P^}  includes  nonminimum-phase 

e 

elements  then  exists  only  if  the  performance  specifications  are  com- 

patible with  the  now  limited  bandwidth  [18]  of  . 

The  same  conclusions  apply  to  the  inner  loops.  Again,  right  half- 
plane P.  poles  pose  no  problem,  but  such  zeros  impose  limitations  on 


L. 


iO  • 


VI.  CONCLUSIONS 

For  a class  of  feedback  systems  with  large  uncertainty,  a multiple- 
loop  design  results  in  sensor  noise  sensitivity  much  smaller  than  in  a 
single-loop  design  satisfying  the  same  specifications.  The  designer  can 
divide  up  the  feedback  burden  among  the  loops  in  a sensibly  optimum 
manner,  wherein  the  uncertainties  of  the  plant  sections,  their  levels  and 
associated  sensor  noise  sources  play  important  roles.  An  important  feature 
of  the  design  techniques  is  its  transparency.  In  return  for  learning  the 
mechanics  of  sensitivity  reduction  in  the  language  of  frequency  response, 
there  is  gained  excellent  insight  into  the  trade-offs  between  the  loops 
and  the  overall  cost  of  design  in  terms  of  bandwidth  and  noise  sensitivity  - 
even  without  performing  the  detailed  design. 
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It  is  discouraging  that  we  must  at  this  time  separately  develop  a 
design  technique  for  each  different  structure.  However,  it  is  encouraging 
that  although  the  present  derivation  is  much  more  difficult  than  for  the 
cascade  system,  the  results  are  remarkably  similar.  This  leads  to  the 
expectation  of  similar  results  for  any  multiple-loop  single  input-output 
structure.  It  is  probably  necessary  to  extend  quantitative  design  to  some 
additional  complex  structures  before  the  general  pattern  will  become  clear 
for  any  multiple  loop,  single  input-output  plant. 
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APPENDIX  I - BOUNDS  (tu)  ON  FIRST  INNER  LOOP 

Sec.  Ill  presented  without  explanation  the  bounds  Bj  (to)  in  terms  of  the 
Intervals  1^,  1^,  1^,  of  L^Cju).  The  explanation  is  available  by  considering 
the  uncertainty  or  variation  set  (9c) 


JvAjJvl 


7/°  - of 


(5a). 


1/^  is  the  line  AB  in  Figs.  A14-(^whereas 

is  the  larger  set  ABC.D.  , a function  of  L,„  and  X.  . The  point  A is 

j j 10  Imax  ^ 

always  the  nominal  L.  , X,  = 1,  k.  = a. , k = a irrespective  of  the  value 

0 1 b b c c 

of  Ljq  , because  that  is  the  objective  of  the  Bj  (oj)  . Attention  is  focused  on 
the  range  -n  < Arg  Ljq  < 0 . The  following  properties  of  are  important. 

(PI)  In  Fig.  Ala,  as  1Ljq|  is  increased  at  fixed  Arg  , boundaries 
BC.  , BD.  shift  dov'jnward  - compare  BC^C^'  at  0 db  with  ^^2^2* 

BC.C,'  at  “AO  db:  and  similarly  the  BD.D.'. 

11  II 

(P2)  For  fixed  L,.  , the  effect  of  increase  in  X,  is  extension  of  the  BC., 

10  I I ’ 

AD.,  i.e.,  widening  of  the  regions  by  decreasing  amounts,  to  a maximum  of 


|0  - tan 


msinB 

l-HncosG 


at  Xj  = <»,  where  = m This  effect  of  large  Xj  is  important  in  explaining 


the  nature  of  B2 (u)  ■ 


(P3)  For  given  Xj^^^  and  |Ljq[  , at  /■\Q  “ ® mirror  image 

(about  AB)  of  ^ ' at  ~ ' 

The  upper  bounds  of  Bj  (m)  in  1^  are  explained  by  property  (Pl)  in 
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Fig.  Ala.  A family  of  at  fixed  ' “90°  is  tried  at  to  = ^lO  e 1^ 

i.e.,  point  A of  |/*  is  set  at  LQ(jAO)  =-32  db  y^l30°  (from  Fig.  3) 
at  which  Fig.  1 requires 


A 


T 

max 


T . 

mm 


< 3^.3  db  . 


At  to  = *»0  (1),  (2)  are  precisely  satisfied.  It  is  seen  in  Fig.  Ala  that  at 

Arg  Ljq  = -90°  , KjqI  < "20  db  is  OK  while  KjqI  > 0 db  is  not  because 

a|t|  = |-3^"2.3|  = 36.3  db  and  larger  |L^q|  gives  larger  aIt|  . The  upper 

bound  here  is  between  0 db  and  -20  db.  From  a study  of  the  shape  of  constant 

|L/1+L|  loci  on  the  Nichols  chart,  it  is  seen  that  this  result  applies  for  all 
toe  1^  at  which  Arg  < -90°  . In  Fig.  3.  there  is  a small  interval  in  which 
Arg  Lq  > -90°  and  in  general  there  may  be  a low  frequency  region  where 
Arg  Lq  > -90°  . However,  the  final  result  is  basically  the  same,  because  of 
the  very  small  sensitivity  of  the  loci  of  constant  |L/1+L]  on  the  Nichols 
chart  at  large  |l|  . 

It  is  worth  noting  that  if  did  not  exist  at  all,  then  V''  i = 1, 

...n  would  only  be  a much  longer  vertical  line  with  lowest  point  at  A . 

From  Fig.  Ala,  both  (1)  and  (2)  would  still  be  satisfied.  Thus  for  (oe  I.  , 

M 

designed  for  P , uncertainty  only,  automatically  handles  P,...P 
u D u In 

uncertainty  as  well.  However  L,  is  needed  in  l„  , precluding  L,  5 0 
In  1^  and  giving  there  upper  bounds  as  in  Fig.  Ala.  Similarly  note  that  in 
1^,  Bj  (oj)  are  hardly  affected  by  large  Increase  of  - see  Fig.  Ala. 

Therefore  Ljq  could  handle  the  entire  uncertainty  of  P^  i.e.  in  place  of 

X,  if  G.J  = G~  = G =0. 

123  n 


A 


2k 


Property  (Pl)  also  explains  in  Fig.  Alb  the  lower  bounds  in  Ij,. 

p 


penetrates  into  M < 2.3  db,  violating  (2).  Thus  i 


m 


Fig.  Alb,  at  Arg  Lj^  = 0,  l*-^olmin  ~ range,  (2)  easily  dominates 

so  there  is  no  danger  of  violating  (1)  (cf  Fig.  Ala)  except  possibly  at 
very  large  Ihol*  which  would  not  be  used  anyhow.  Here  too,  Aj  could  be 
increased  to  " without  affecting  Bj((jj)  - recall  (P2),  the  effect  of  large 


Aj  on  in  Fig.  Ala  and  the  critical  factors  in  Fig.  Alb.  Thus  there  is 


no  need  for  in  Ig  as  well.  (PI)  also  explains  in  Fig.  Ale  the 

upper  and  lower  bounds  in  Ig.  Thus,  at  = 6.  either  < m2 

or  > mr.  From  (P2)  the  width  of  is  < [Arg  Lj^l-  Hence,  Fig.  Ale 


>i^«  11V./111  \i^/  ciici  rfivjcii  wi  ^10 

shows  that  as  to  increases  in  1^.,  the  value  of  -Arg  •■lo  for  which  al 1 11..]  are 


acceptable,  increases  steadily,  explaining  why  the  B,  (to)  closed  curve?  shrink 


to  Bj^  in  Fig.  6b.  Bj^  length  is 


1 


^1 


db 


because  at  Arg  L 


10 


IT  (say  at 


3] 

1,000  here)  1 + AjLjq  = 1 - AJLjqI  with  |Ljg|  < l/Aj^^^g^  = -g- ^ otherwi se 


'110 


'1>"10' 


Imax 


.1 


y'  extends  in  length  to  «»  and  being  360“  wide,  must  intersect  with  the 
forbidden  |j^l  < Y “ 2.3  db  regions  located  at  Arg  L = ± nir,  n = 1,3,... 


This  is  also  seen  from  (10),  for  let  L.  = A.L.  , and  P = P..  for  a i . 

' I 110  a ao 


Then  at  Arg  L.  =-Tr,  |L.  |<1/A.  is  essential , otherwise  Ic./D. I is 
^10  ’ ' 10'  imax  '.11' 

inf  ini te  at  A. 

imax 

Increase  of  Aj  affects  the  bounds  at  m2,  requiring  iI-jqI  < 3m2, 

3 < 1,  but  not  the  lower  boundary  at  m^.  In  Fig.  6b  it  is  seen  that  Ljg 
lies  on  the  upper  part  of  Bj (w)  for  most  of  Ij^,  so  Lj^  designed  to 


handle  Pj  only,  can  also  cope  with  f’2’"''^n  L2=...=L^=0. 


il' 


if 
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APPENDIX  2 - BOUNDS  ON  SECOND  AND  HIGHER  INNER  LOOPS 

The  function  of  is  to  guarantee  that  Ljq  is  satisfactory 

despite  its  design  on  the  basis  of  P^^  of  (9a).  It  was  seen  in 

Appendix  I that  for  w e Ig  snd  part  of  ■c’  ho  suffices  i.e., 

L^q  may  be  zero.  This  is  so  only  for  co  < at  which  ~ 

It  was  noted  also  that  I I 1 — ' *"2  needed  for  oj  = 

''imax 

Hence,  L^q  =0  is  impossible  in  Ij^  and  it  is  not  surprising  that 

the  B2(w)  there  are  upper  bounds  (recall  in  Appendix  1 precisely 

the  same  situation  for  LjQ(ja))  in  I^).  At  oj  = Wj^eljg, 

£>0  (.38  in  the  example),  so  the  denominator  (12a)  of  P^  is 

^30  it  0 at  Arg  L2q=0,  |L2oI>  

2 max 

= 1.^  here.  So  there  is  a lower  bound  on  I *-20 ^j^lTr  ^ I which  is  a function 
of  Arg  L2Q. 

2 

To  find  82(0))  in  Ijj.  it  is  necessary  to  use  P^  of  (12a)  in 

1 2 
place  of  P^  of  (9a).  It  is  convenient,  however,  to  express  P^  in 

1*  1 

terms  of  P^  , defined  as  with  A|^2  replacing  Xj,  because 


^2®^20’ 


QC|  QD.  * O.C.  . + aXjLjd  I * ^ 

Wi  " 1 - -QD,  + D,V  - at  aL,,„  “ I + 4,, 

' * ^24o 

OC.  = oa  + QCj  = >‘,>^2'-|0  ■"  -T+~L~q  ' ’ 


Fig.  A2  was  sketched  for  u = 2000,  Xj  = kO,  which  (Fig.  6b) 

LjO  “ .0158^/230°,  = -127  db  y^h30° , for  assumed  Arg  L^q  = -117°  constant. 

The  D.  describe  an  arc  of  a circle  as  1L2q|  's  varied,  as  do  the  C.  drawn  for 
= 10  = max,  i.e.  DjCj  = 10  D.V.  Clearly  for  K^qI  < < I,  OV/OC  ->  1 
and  for  iL^p]  > > 1,  lOV/OCl  < 1,  so  such  |L2q|  are  acceptable.  Obviously 
3 < A2^3^»  5 resulting  C.  circle  passes  through  0,  giving  infinite  OV/OC. 

and  the  resulting  ^2/^{P^}  passes  thru  M = 2.3  db.  Thus,  3 upper  and  lower 
bounds  in  this  oj  range.  As  w increases,  KjqI  and  its  angle  drecrease,  so 
the  arc  C.Cj...  does  not  extend  to  0 in  Fig.  A2  and  any  11-2q|  is  acceptable. 
Hence,  the  82(0))  tend  to  a line  at  -tt,  from  0 to  (32^*^2^db" 
shown  in  Fig.  7.  including  a sensibly  optimum  L2Q(ju))  with  its  intervals 

'2A’  ‘28’  '2c' 

For  the  third  inner  loop  (if  n > 2),  is  needed  and  there  is  an 
analogous  situation  with  respect  to  1^^.  At  l2g>  ^277  which  Arg  = - it) 
is  very  large  (''^  65OO)  and  as  before,  there  is  a lower  bound  on  L^q  at 

‘^2tt*  “-“2ti’  l‘■lol  so  (in  13d),  Denom.  (P^)  1 + 

+ AjA21-2q(1  + 1-jq)  » similar  to  Denom.  (P^)  if  i is  replaced  by  i - I . 

3 2* 

P^/P^  similar  to  (Al)  is  obtained  giving  a figure  similar  to  Fig.  A2  and 
analogous  8^(0)).  The  process  is  continued  to  B^(u)),  ...,  B^(a)). 
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20  AB^y^  ACT  (Continue  on  reverse  side  it  necessary  and  identify  hy  black  number) 

There  is  given  a single  input-output  linear,  time-invariant 
plant  with  large  parameter  uncertainty  consisting  of  two  parallel 
branches,  one  of  which  has  n internal  sensing  points.  The 
objective  is  to  satisfy  specified  frequency  domain  bounds  on  tho^ 
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system  response  to  commands  and  disturbances  over  the  parameter 
range,  and  to  do  so  with  sensibly  minimum  net  effect  at  the  plant 
input,  of  the  n + 1 sensor  noise  sources.  The  basic  problem  is 
how  to  best  divide  tlie  feedback  bvirden  among  the  n T 1 available 
feedback  loops  lij^.  The  procedure  developed  lias  high  transparency, 
giving  early  perspective  on  the  loop  bandwidths,  permitting 
approximate  loop  trade-offs  without  a detailed  design,  n Vilhile 
[the  development  is  more  difficult  than  in  the  single  cascaded 
plant  system,  the  procedure  and  final  results  are  very  similar: 
Each  has  only  one  distinct  frequency  range  say  in  which 

there  is  trade-off  between  L.  and  L.^.,  and  w • . i with 

steadily  increasing  loop  bandwidths  going  backwards  from  plant 

I 

output  to  input. It  is  shown  that  for  a class  of  problems  the 
sensor  ooise  effects  can  be  tremendously  reduced,  when  compared 
to  an  optimum  single-loop  design  satisfying  the  same  specifica- 
tions. 
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